Control of bound-pair transport by periodic driving 
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We investigate the effect of periodic driving by an external field on systems with attractive pair- 
ing interactions. These include spin systems (like the ferromagnetic XXZ model) as well as ul- 
tracold fermionic atoms described by the attractive Hubbard model. We show that a well-known 
phenomenon seen in periodically driven systems — the renormalization of the exchange coupling 
strength — acts selectively on bound-pairs of spins/atoms, relative to magnon/bare atom states. 
Thus one can control the direction and speed of transport of bound-pair relative to magnon/unpaired 
atom states, and thus coherently achieve spatial separation of these components. Applications to 
recent experiments on transport with fermionic atoms in optical lattices which consist of mixtures 
of bound-pairs and bare atoms are discussed. 
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I. INTRODUCTION 

The dynamics of quantum spin transport and in partic- 
ular spin correlations are of much current interest in the 
field of quantum information [1]. Cold-atoms in optical 
lattices provide clean realizations of a range of many- 
body Hamiltonians, including well-known spin systems. 
There is also enormous interest in pairing phenomena, 
motivated by many ground-breaking experiments with 
ultracold fermionic atoms in optical lattices Q. 

Neglecting particle interactions, many simple many- 
body Hamiltonians, such as Heisenberg, Hubbard or 
Bose-Hubbard, comprise a hopping term —JH^ charac- 
terized by an exchange coupling strength J. If such a 
system is subjected to an additional spatially-linear, but 
oscillating field, we have a total time-dependent Hamil- 
tonian: 

N 

H(t) = -JH h + B8mut^2na*/2. (1) 

n=l 

For these, as well as other analogous driven systems, 
there are regimes where the exchange strength takes an 
effective, renormalized, value 0, 0, [B|, : 

Jeff = JJb(-), (2) 
UJ 

where Jq denotes an ordinary Bessel function. The inter- 
site transport is completely suppressed at values of ^ = 
2.4, 5.52, ... corresponding to Jo(^j) = 0. The oscillating 
potential in Eq.([T]) can be implemented with ultracold 
atoms in shaken optical lattices; thus Eq.([2|) was recently 
demonstrated experimentally [7|, [8[ by mapping | J e ff | as 
a function of the ratio — . 

UJ 

In the high frequency regime, where uj ^> J, the sup- 
pression is often termed Coherent Destruction of Tunnel- 
ing (CDT) 0, The behavior of the underlying quan- 
tum spectrum and the application of CDT to controlling 



the superfluid-Mott Insulator transition in atomic sys- 
tems were investigated theoretically in Ref. [9|, |_10(] . 

However, this CDT is also closely related to a phe- 
nomenon termed Dynamic Localization (DL) — identified 
even earlier in 1986 [3] — which also impedes transport 
in periodically driven systems at »Jo(fj) — 0. It too has 
been theoretically investigated in atomic systems [ll|, [l2| . 
The precise relationship between CDT and DL remains 
of much interest: a recent theoretical analysis is provided 
by Ref. [l3j], but in brief, the CDT mechanism is associ- 
ated with high-frequency driving (uj ^> J) and complete 
suppression of hopping. For CDT, particles remain com- 
pletely frozen at their original sites even for a 2-site sys- 
tem: CDT was initially identified in driven double-well 
systems 0, . CDT can persist even in the presence of 
some inter-particle interactions [9|, DL, on the other 
hand, entails a less complete suppression of transport: 
for Jo(fj) — 0, the single-particle wavepacket position 
may oscillate, but the particle returns periodically to its 
original position. DL is associated with lower frequency 
driving, negligible particle interactions and the large TV 
limit jjj. 

Here we investigate for the first time periodic-driving 
for systems with a bound-pair component. In the absence 
of driving, transport for fermionic systems with attrac- 
tive pairing interactions has been previously investigated 
in, e.g., ferromagnetic XXZ spin models [14| and even ex- 
perimentally with ultracold atoms which are a mixture 
of bound-pairs and bare atoms, in a regime of the at- 
tractive Hubbard model [IBj]. Our main finding is that 
the DL regime (but not CDT) can provide a mechanism 
to spatially separate the paired and unpaired fraction in 
such spin or cold atoms systems. This is because the 
driving can generate a hitherto unnoticed mechanism of 
directed motion; in contrast, previous studies to date, 
such as Refs. 0? S EDI ? probed only |J e ff|, i.e., the pres- 
ence or absence of diffusive expansion of an atomic cloud, 
not global transport. 
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We note that both the DL and CDT transport- 
suppression mechanisms are usually analyzed in terms 
of the stationary states (Floquet states) of the driven 
Hamiltonian: typically, the Floquet eigenstates (i.e., 
their quasienergies) become degenerate or approximately 
degenerate if Jo(^j) — 0. However, here, in the N — > oo 
limit, we analyze the quantum transport without any ref- 
erence to the detailed quantum structure, or quasienergy 
near-degeneracies: we obtain the renormalization Eq.([2|) 
(as well as the ratchet mechanism) from a straightfor- 
ward integration of the classical equations of motion of 
the driven Hamiltonian. This entails a somewhat mod- 
ified perspective, since now the key division is not be- 
tween high-frequency CDT freezing-of-motion and low- 
frequency DL wavepacket revival. It is rather between 
the N — » oo limit, where both high and low frequency 
suppression of transport arises as a classical phenomenon 
(i.e., it is purely an effect in the group velocity), and 
the small N limit, where a quant al analysis of a few 
Floquet state degeneracies remains essential — and where 
only high frequency driving suppresses transport. The ef- 
fect of inter-particle interactions seems difficult to treat 
classically. But, below we identify an example where we 
can include the effect of certain important interactions 
by considering two limiting effective Hamiltonians. 

The paper is organized as follows: In Sec. HI we intro- 
duce our models. One is a Heisenberg XXZ spin chain; its 
eigenstates include both magnon-like spin-waves (analo- 
gous to Bloch waves) and bound-pair states. Explicit 
expressions for both are obtainable via the Bethe ansatz. 
The close correspondence with the two-particle attractive 
Hubbard model, which also includes both bare atom and 
bound-pair states, is explained. The Hubbard model is 
of particular interest as it is already realized in current 
cold atoms experiments. In Sec. IIII( the quantum CDT 
mechanism, dominant for low N and co ^> J is intro- 
duced. Its effects are demonstrated for both unpaired 
and bound-pair states by exact numerical solution of the 
quantum Hamiltonian. In Sec. HV( the low- frequency, 
high- TV regime is investigated; control of wavepacket dy- 
namics is shown for the one and two-particle regimes. 
The underlying mechanism of directed motion are ob- 
tained from the classical equations of motion. The be- 
havior of both bound and unpaired cases are analyzed 
with two limiting classical Hamiltonians and are shown 
to agree with the full quantum solutions. Conclusions 
and outlook are given in Sec. [VI 



II. STATIC MODELS 

In general, inter-particle interactions are important so 
the full many-body Hamiltonian takes the form JH = 
JH n + (J A) H- m t, where H[ n t is an interaction term char- 
acterized by an interaction strength U = J A. For in- 
stance, the well-known spin-1/2 Heisenberg XXZ ferro- 
magnetic chain of length N, is governed by the Hamilto- 



nian: 

N j JA 

~JH= -^[^{°+(T n+1 +(T n CT+ +1 ) + — (3) 
n=l 

where n G [1 : N] indicates the n-th site, and A denotes 
the anisotropy. For large A (and low density of excita- 
tions), there are two dominant classes of quasi-particle 
states: magnon-like and bound-pair states. 

The Hamiltonian Eq. ([3]) conserves the number of spin- 
flips; a single excitation represents a spin- wave, or 
magnon, which distributes a single spin-flip through- 
out the chain. The spin- wave eigenstates of Eq.((3]) are 
| ft) = -±= J2^=i e%nK \ n ) (f° r P er i°dic boundary condi- 
tions) where \n) denotes a state with the spin- flip at sites 
n; they obey the dispersion relation E K — Eq = — J cos ft, 
where Eq is the ground state. Single magnon transport 
is thus analogous to free atoms (Bloch waves) in the low- 
est band of a lattice. Higher excited states correspond 
to multiple spin-waves which interact when they coin- 
cide through (i) an exclusion process (no two spin flips 
can simultaneously occupy the same site) and (ii) an ef- 
fective attractive interaction induced by the (JA)H[ n t = 
E n T a n a n+i interaction term. For the long wavelength 
processes considered in our key results Figsl2]and[3l only 
(ii) is of significance. 

Now we consider the two-particle case. For the two- 
excitation case, the eigenstates of Eq.(j3]) may be ex- 
pressed, via the Bethe ansatz, as spin waves ful. fl6j|: 

\ki,*2) = A KuK2 a niZ \ni,n 2 ), (4) 

n 1 <n 2 <N 

where A KliK2 is a normalization constant and 1 721,712) in- 
dicates a spin-flip at sites n\ and Further details 
are in Refs. [3, G3] but the eigenstates divide into 
two distinct classes: (1) Magnon-like scattering states, 
where the spins move separately, with dispersion rela- 
tion E KljK2 — Eq = J (2A — cosfti — cosft2) which is a 
straightforward extension of the singly-excited case and 
(2) Bound-pair states, for which the probability ampli- 
tudes decay exponentially with the separation of the flips. 
For A > 0, and as N — > 00 they obey the dispersion rela- 
tion, E Kl , K2 — Eq = JA — ^ [1 + cos(fti + K2)] • Note that 
the sum (fti + ^2) is always real. In Ref. [14( a perturba- 
tion expansion in spin coupling strength J is employed to 
produce an effective Hamiltonian when A ^> 1. In this 
approximation the bound state amplitudes are 

CL n 1 ,n 2 ~ °ni,n 2 -l e \ D ) 

and thus the bound states become confined to the near- 
est neighbor (NN) subspace, {|n, n+1)}; their dispersion 
relation remains unchanged from the above. In "center 
of mass" coordinates 2ft = fti + ^2, the bound states are 
given by |ft) = A K J^^Li e 2inKj |n, n + 1) and are of anal- 
ogous form to the single magnon solution. Namely, two 
initially neighboring spin-flips thus hop together but their 
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FIG. 1: (color online) CDT demonstrated in a XXZ chain 
with N — 20 and large A = 8. For A ^> 1, initial spin- states 
of adjacent spin flips \ip(t = 0)) = |n, n+ 1) overlap only with 
bound-pair states; initial states with well- separated spin-flips 
\ip(t = 0)) = |ni,7i2) = |5, 15) overlap only with magnon-like 
scattering states. The fidelity F(t) = \(ip(t = 0)\ip(t))\ 2 is 
shown as a function of time. Insets plot T90, the time taken 
for the fidelity to fall to 0.9, so show the CDT resonances, (a) 
J' — 1/8. CDT for magnon states. Spins freeze at zeros of 
Jo(B') (B f ~ 2.4,5.5...) but delocalize elsewhere, (b) ^ = 
1/8. CDT for bound-pairs. Adjacent spins freeze at fields for 
which J (2B') — (B' c± 1.2, 2.7...) 



propagation speed is slower relative to a single spin- flip: 
transport occurs by a second order effective Hamiltonian 
and an effective coupling J e ff = ^ < J. 

The transport is rather analogous to that seen in the 
attractive Hubbard model —JH = J a (Cka c (k+i)a + 
H.c.) + J A n^riki characterized by tunneling ampli- 
tude J and interaction U = JA; for U < 0, local bound- 
pairs form which tunnel by a second-order process, with 
amplitude of order ~ J 2 /U ~ J/ A akin to the spin sys- 
tem. This situation was investigated experimentally us- 
ing fermionic atoms for U < in Ref. [15[ ; note of course 
that the atomic local bound-pairs (BPs) correspond to a 
pair of atoms of opposite spin occupying the same lattice 
site, while the spin BPs occupy adjacent sites. 



III. QUANTUM DYNAMICS: HIGH 
FREQUENCY DRIVING AND COHERENT 
DESTRUCTION OF TUNNELING 

We now consider the regime with B,co ^> J which can 
lead to coherent destruction of tunneling. Using units 
of scaled time t' = cot, the quantum dynamics generated 
by Eq.(pQ) depends only on two global parameters J' = 
J/co and B' = B/co (for a given A). The high-frequency 
driving condition for CDT thus becomes J' « 1. We 
follow the usual procedure for justifying CDT [H, [To| : 
for B' ^> J', we initially neglect and set H(t) = 
S'sini'£Li«</2- 

Then, the Schrodinger equation has the following time- 
dependent solutions: 



|ni,n 2 ,t) 



= p [iB' (m +ri2) cos t] 



|ni,n 2 ), 



(6) 



to within a constant phase term e^ B ^ cos£ ] where up 



^2n=i n = N(N + l)/4. The stationary states of peri- 
odically driven systems (Flo quet states) have a Brillouin 
zone type structure [9, 10] and the above correspond to 
a family of solutions |ni,ri2,m, t) = |ni, 77,2, t)e irnt which 
are periodic in time. The multi-photon band index m 
is used to label states which differ only by the absorp- 
tion or emission of m quanta of energy E m = mhco. The 
next step is to treat the spin-exchange as a perturbation, 
in the extended Hilbert space where time appears as an 
extra coordinate. The matrix elements of are simply 

((n[ 1 n 2l m f 1 t\JH h \n ll n 2l m 1 t}} ~ £ mm ' 
J Jo - m) + B'(n' 2 - n 2 )} (n' 1; r^\H h \rn, n 2 ). (7) 

Here ((.|.)) = ^ f Q T (-\-) denotes a scalar product 
in the extended space [9] and (.|.) denotes the scalar 
product in position coordinates. It is easily seen that 
in general, the matrix elements (n[, n 2 |i?h|^i, n 2 ) = 

\H h \n 1 ,n 2 ). In 

other words, they involve hopping of single spins only. 
Thus, we have 

J J [B , (n , 1 - ni ) + B\n' 2 - n 2 )\ = JJ {B') (8) 

as in Eq.Q. This also corresponds to the situation ex- 
plored in the cold atoms experiments @, Q. A more 
careful analysis [9[ reveals terms coupling off-diagonal in 
the photon band index m. If these off-diagonal terms are 
significant, the above analysis is no longer valid; however, 
it has been shown that provided the energy scale for the 
band separation far exceeds the intra-band energy scale 
(i.e., co ^> J), these couplings are negligible. Thus within 
the above framework, the Eq.([2j) renormalization holds 
for high-frequency driving. 

However, the bound pairs, for A ^> 1, evolve under an 
effective second-order Hamiltonian H e ^ [lj] which results 
in the two spins hopping together. When n 2 = n\ + 1, for 
instance, the argument of the Bessel function becomes 

B'K - n i) + B '[ n i + 1 - K + !)] = 25 'K - n i)- ( 9 ) 

Thus, the transport is then determined by matrix ele- 
ments ^ K Jo(2B'){n,n + l\H e ^\n ± l,n + 1 ± 1): note 
the doubling of the argument of the Bessel function. 
Since the Bessel is an oscillating function one can then, 
for example, choose a value of B' for which J e R 
JJo(B') for the magnon-like states is positive, while 
J eff ~ t^Jo(2B') for the bound-pairs is negative (or 
zero). 

An arbitrary spin state, in general, has a projection 
on both the magnon-like and the bound-pair eigenstates. 
For A ^> 1, however, one may prepare a good approx- 
imation to a pure bound-pair state by simply flipping 
two adjacent spins of a ferromagnetic chain in its ground 
state (all spins aligned). Conversely, two well-separated 
spins will approximate pure magnon-like states. In a 
corresponding atomic experiment with Fermionic atoms 
one may either consider unpaired single atoms or a pure 
bound-pair, as well as a superposition of these two ex- 
tremes. 
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In FigJTJwe demonstrate CDT for the two extremes (an 
initial state which is either a pure magnon state or a pure 
bound-pair for large A). Driving with a field B' = 5.52 
"freezes" two initially well-separated spin-flips at their 
original sites, since J${B f ) = 0. In the absence of driving, 
or at values of B' with J'o(B') ^ 0, both spins rapidly 
diffuse along the chain. On the other hand, if the two- 
flips are initially adjacent, they remain frozen at their 
positions if Jq(2B') = and delocalize otherwise. 

Other than for these extremes, CDT (e.g., for initial 
states which are superpositions of magnon/bound-pairs) 
is less effective; while both the magnon-like scattering 
states and bound-pair states are eigenstates of the Hamil- 
tonian Eq. (|3]), the respective subspaces are coupled by 
the driving: in general, for CDT, the driving strength 
B ^> JA, where J (A — 2) is the energetic separation 
between the two subspaces. 



IV. THE LARGE- REGIME AND DYNAMIC 
LOCALIZATION 

We now investigate the regime N — > oo, usually as- 
sociated with Dynamic Localization. We can work in 
the more favorable regime JA > B where coupling be- 
tween bound-pairs and unpaired states is suppressed. We 
evolve Eq.(pQ) for an open chain with N = 100 sites (for 
computational reasons, much larger N is difficult). Since 
our analysis below does not introduce any high-frequency 
condition, we can also consider low-frequency driving, for 
which B',J' > 1. 

We prepare a one- spin flip Gaussian wavepacket, at 
the center of the chain. Figure [SJa) shows that the spin- 
packet unexpectedly moves as a whole, and with little 
spreading, along the chain, but its direction of travel re- 
verses at B' c± 5.5: for B' = 5.3 it moves upwards, for 
B' — 5.7 it moves downwards; for B' = 5.5 it exhibits 
behavior characteristic of Dynamic Localization (static, 
but with oscillations). The results for B' = 6.3 may seem 
at first sight even more puzzling: the wavepacket' s cen- 
ter of mass is static, but the excitation diffuses along the 
chain. In order to understand these results, one notes 
that Fig|2](a) demonstrates a combination of DL type 
renormalization, as well as a (hitherto unnoticed) type 
of ratchet (meaning directed motion without bias). 

We note that N plays the role of an effective H. Thus, 
in the presence of driving and in the TV — > oo limit, one 
may map the system to an "image" classical Hamiltonian: 



H(x,p) = — J' cosp — B'x sint' 



(10) 



by mapping site to continuous position n — > x and n — > p. 
Integrating Hamilton's classical equations of motion over 
one period, it is easy to see that the distance traveled per 
period v = (x(t + T) — x(t)) j (27r), where T = 2tt/co, for 
a particle with position x{t = 0) = xo and momentum 
p(t = 0) = po is simply: 
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v = J'J (B')sm( Po + B'). 



(11) 
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FIG. 2: (color online) (a) Probability distribution |^(n,£)| 2 
for an initially Gaussian wavepacket containing a single exci- 
tation, near the dynamic localization (DL) regime, for J' — 8. 
Two zero-velocity cases are shown: B' = 5.5 corresponds to 
actual DL, with Jo(B') = and suppression of hopping; while 
B' — 6.3 corresponds to sini^ = 0, so suppression of directed 
motion, but not of hopping: the excitation diffuses along the 
chain, (b) Average velocity of the wavepackets as a function 
of B' . Numerics were obtained from plots like those in (a). 
They are well- fitted by the magnon velocity Eq. ([ll[) , while the 
predicted bound-pair behavior corresponds instead to Eq. (fT2|) 
shown here for A = 2. The vertical lines a,b,c correspond to 
the multiply excited packets shown in Fig. [3] 



In effect, this is the center of mass velocity of the 
wavepacket. It is independent of initial position x$. The 
renormalization of the velocity already appears, but it 
is no longer a question of simply renormalizing J' by 
Jq(B') because there is also a sin(po + B') factor which 
can also change sign or suppress transport. We find, how- 
ever, that it does provide great advantages in reducing 
dispersive spreading. 

An initial Gaussian wavepacket, initially localized over 
a finite number of sites 1/S n , corresponds to local- 
ization in momentum of ~ 5 ni in the image phases- 
pace. Thus the corresponding momentum distribution is 
N(p) ~ exp[— (po — (p)) 2 /^n}' Low-energy spin wavepack- 
ets (and cold atom clouds) correspond to distributions 
well- localized about zero momentum, i.e., n ~ thus 
(po) ~ 0. The distribution thus samples velocities ~ 
J'J (B')sm(B'±5 n ). Hence provided mir < \B'±5 n \ < 
(m + 1)tt, all momenta —S n <po<S n correspond to 
the same direction of motion and there is little disper- 
sion. This results in directed motion. In contrast, when 
B' c± mir with m = 1,2,..., trajectories with po < 
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Time (t/T) n 2 

FIG. 3: (color online) Dynamics of spin states which are su- 
perpositions of bound-pair and magnon states, for param- 
eters shown in FigEfb). (a): J' = 10, A = 2,B' = 2.4. 
The magnon states are stopped, while bound-pairs travel to 
the end of the spin-chain, (b) J' = 10, A = 2,B' = 2.53. 
Bound-pair and magnon states move in opposite direction 
(bound-pairs move upwards), (c) J' — 10, A = 2, B' — 
5.34. The faster magnon packets are slowed by the field 
so magnon and bound-pair speeds become equalized and 
they now travel together, (d) J' = 2, A = h,B' = 4.33. 
Bound-pairs remain static, magnons travel, (e) Spin-spin 
correlation, | (m, ri2 | ^(t))| 2 , at t/T = 10 for the spin dynamics 
of (d). The bound-pair component remains stationary, in its 
initial position at the center, while the magnon-like compo- 
nents travel to the end of the chain. 

move in opposite direction to those with po > and the 
wavepacket, though static, spreads diffusively along the 
chain. In the presence of DL (B f = 5.5), Jo(B f ) = 0, so 
there is, of course no transport. 

In Fig|2](b) we see that a calculation of the veloc- 
ity from the numerics [displacement of the quantum 
wavepacket per period obtained from a quantum solu- 
tion of Eq.JT}] is in excellent agreement with Eq. ([TT]) . 

Conversely, a single bound-pair state moves under 
a different image Hamiltonian H = cos 2P — 
2B'Xsmt', where 2P = p 1 + p 2 and X = {x x + x 2 )/2, 
indicate center of mass coordinates. Thus the bound-pair 
velocity: 

v hs = ^M2B') sin(p + 2B'). (12) 

For multiple excitations, we analyze the numerics by as- 
suming that the magnon-like projection of the wavepack- 
ets obey Eq. ([TT]h while the Hilbert space fraction of 
bound-pair states moves under Eq. ([T2|) . We thus assume 
that coupling between the two subspaces is negligible, a 
reasonable assumption for B <C J A. 

Figure [2jb) (broken line) shows the corresponding ve- 
locity. One sees that the bound-pair and magnon ve- 
locities may even have opposite signs for the same driv- 
ing field. This enables us to steer the two components 
in opposite directions or to stop one and transport the 
other. These different situations are demonstrated in 
Fig|3l where the initial state is the product of two Gaus- 
sian wavepackets peaked on k ~ and at positions 
m = 45 and ri2 = 50. The k ~ condition corresponds 
to the low-energy distribution — also typical of an ultra- 
cold gas (centered on po ~ 0). There is no special signifi- 
cance in the initial chosen positions [neither Eq. (fTT]) nor 



([T2]) depends on initial x] other than choosing an initial 
state with n\ and n 2 close but not adjacent yields an 
appreciable projection in both bound-pair and magnon 
subspaces. Note also that the extent of the initial distri- 
bution should not exceed that of the chain/lattice, so that 
center of mass displacement is observable. Equivalent ly, 
in a Hubbard Hamiltonian, two delocalized atoms with 
distribution peaked one or two sites apart would yield 
an appreciable probability to both form a pair /remain 
unpaired. 

The directed wavepacket motion is in sharp contrast to 
the usual motion seen in studies of CDT/DL transport; 
also to the undriven case: while magnon states diffuse 
faster than bound-pair components (by a factor of 2A) 
both ultimately simply delocalize along the length of the 
chain. Figure [3je) illustrates the spin correlations for 
FigE9(d) (for large A = 5) where the bound-pair compo- 
nent is immobilized, while the component with separated 
spins is transmitted. 



V. CONCLUSIONS AND OUTLOOK 

In conclusion, we have investigated the transport in 
a system with attractive pairing interactions periodically 
driven by an external field. The effects of periodic driving 
act selectively on unpaired and bound-pair states, which 
enable us to control separately the relative direction and 
speed of each wavepacket of the two states. 

One should also consider higher numbers of particles. 
For low-numbers of excitations (modest filling factors in 
the atomic case) the dominant processes are the above. 
Nevertheless, the effect of including higher numbers of 
particles remains an important question. In the case of 
high-frequency driving and CDT, this is less of a concern: 
e.g., in Ref. [10] it was shown that a 2-particle model 
adequately described the CDT behavior of a full many- 
body Hamiltonian. However, for low frequencies, this 
remains an open question. One of our key findings is the 
onset of center of mass motion in the DL regime, which 
would open new experimental possibilities, and is well 
within reach of current techniques. The experiments of 
Ref. [7] observed the effects of Eq. ([2|) in the range J' ~ 
3 — > 1/30; for lower frequencies, the effect was lost due 
to inter-particle interactions. This range of J' overlaps 
with our DL wave-packet splitting regime [e.g., Fig|2jd) 
corresponds to J' = 2 and U/J = 5]. In addition, a 
system corresponding to the attractive Hubbard model 
is somewhat more favorable than the bosonic atoms used 
in CDT experiments since there are two good DL limits 
provided by the fully paired and unpaired extremes, while 
for the bosonic systems, low-frequency DL occurs only 
in the limit of negligible interactions. A full numerical 
study of the DL regime (high N and many particles) is 
numerically challenging, but future calculations will need 
to address this question more fully. 

Within the regime of validity of the present work, even 
for not too large J A /B we see that the wavepacket splits 
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cleanly into two components, whose relative speed and 
direction may be controlled as in Eqs. (fTT]) and ([T2]) . This 
shows that coupling between the magnon and bound-pair 
subspaces by the external driving is limited. In fact it is 
easy to show that, once the wavepackets separate, fur- 
ther interaction is negligible. If the separation is slow, 
collisions between the paired and unpaired particles may 
be important; but this itself exemplifies a topic of much 
current interest; for example recent experiments probed 
transport of impurity wavepackets accelerating through 
a static cold cloud by means of a linear potential V = gx 
due to gravity [18[. The study of dynamics of bound- 



pairs moving through a cloud due to the linear oscillat- 
ing lattice potential V = Fx cos cot is also within current 
experimental capabilities and offers new possibilities for 
transport in the presence of interactions. 
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